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Abstract. The problem of the flux to a spherical trap in one and three dimensions, 
for diffusing particles undergoing discrete-time jumps with a given radial probability 
distribution, is solved in general, verifying the Smoluchowski-like solution in which the 
effective trap radius is reduced by an amount proportional to the jump length. This 
reduction in the effective trap radius corresponds to the Milne extrapolation length. 
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p{r,t) = — 



for r > i? , (1) 



1. Introduction 

A classical problem in diffusion and reaction theory is the flux of particles, initially 
uniform in space, to a trap. When the particles are moving by simple diffusion, this 
flux can be found by solving the diffusion equation with an adsorbing boundary at the 
surface of the trap. For example, for a system with a spherical trap of radius i?, the 
time-dependent density p(r, t) evolves via the diffusion equation, dtp = DV^p {D being 
the diffusion coefficient), starting from the initial condition p(r, 0) = po for r > R and 
p(r, 0) = for r < R. The diffusion equation is subject to the boundary conditions, 
p{r,t) = for r = i? and p{r,t) — >■ po as r — oo, for all time t. For spatial dimensions 
d > 2, the density profile outside the sphere approaches a steady state as t oo, 
reflecting the fact that Brownian motion is transient for c? > 2 pp. For example, in three 
dimensions, the full time-dependent density profile (which is spherically symmetric) can 
be explicitly obtained for all t, 

i?erfcf^y 

where erfc(2;) = (2/yvr) e~^^ du. In particular, as t — oo, the density approaches a 
stationary profile outside the sphere, 

p(r,oo) = ^(r-i?) . (2) 

The flux $(t) = 4'KR'^D[dp/dr)\r=R, defined as the number of particles falling into the 
trap per unit time, can be easily obtained from p(r, t) in Eq. (P) and one gets the classical 
result [aEl, 

R 

On the other hand, when the particles are moving by discrete jumps at every time 
step r (as opposed to undergoing continuous-time Brownian motion), the above results 
for the stationary density profile in Eq. (0) as well as that for the flux in Eq. (jS)) become 
modified. In Ref. jlj this discrete-jump problem was studied by an iterative numerical 
scheme, for particles undergoing fixed-length jumps of length i and random angles (the 
so-called Rayleigh flights). It was found that in d = 3, while the density profile Pnir) 
after n time steps approaches a stationary limit as — > oo, the stationary profile Poo{f) 
is different from its continuous-time counterpart in Eq. (0). Very far from the surface 
of the sphere, the stationary profile behaves as 

Poo(r) ^ ^ [r - (i? - e)] for r > i? , (4) 

where e = c^ and the dimensionless constant c = 0.29795219... was determined 
numerically On the other hand, on the surface of the sphere, the stationary density 
profile approaches a positive value 



<^{t) = 4:71 RDpo 



1 + 



(3) 



Poo (i?) ^ 0.408245^ , (5) 
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where the constant 0.408245 was evaluated numerically in Ref. jlj. This is in marked 
contrast to the continuous-time Brownian case where the stationary density profile 
vanishes on the surface of the sphere. The distance e = ci is the 'Milne extrapolation 
length' m that represents the distance inside the surface where the far steady-state 
solution in Eq. (j^ extrapolates to zero. The same length e also appears in the expression 
for flux H] 

R-e 



$(t) ^ AnDpoiR-e) 



1 + + o{t-'/') 



(6) 



for large t. This result was obtained numerically in Ref. |3]. 

Recently we have rigorously proven the result in Eq. for the large-distance 
stationary profile and found an analytical expression for c [8 : 



1 [°° dk ^ 
c= — I p In 



TT JO 




, , 0.29795219028... . (7) 
k ' 

Remarkably, the same constant also appeared in the leading finite-size correction to 
the expected maximum of a discrete-time random walker moving on a continuous line 
where at each time step the walker jumps by a distance ^ that is chosen from a uniform 
distribution OE]. In Ref. the constant was evaluated numerically by summing a 
rather complicated double infinite series. Later an exact expression as in Eq. ((Zj) was 
derived in Ref. jTUI. In Ref. [S], we showed why the same constant also appeared in 
the three-dimensional discrete flux problem. For the 3-d flux problem, it was further 
proven jS] that the density on the surface of the sphere approaches a constant as in Eq. 

Poo(r) = -^-. (8) 

Thus the constant 0.408245 in Eq. (0), found numerically in Ref. was proven to be 
simply 1 / -\/6. Indeed, in Ref. jH] , an analytical expression for the full stationary density 
profile for all r > R (or rather its Laplace transform) was obtained. 

Having obtained the stationary density profile, the next objective is to prove 
the expression for the time-dependent flux in Eq. (jH)) that was found numerically in 
Ref. 13]. If one uses the form in Eq. Q of the steady-state solution far from the 
sphere in the expression for the flux in the continuous-time diffusion problem, one 
gets |H| $ = Anr'^D {dpoo{r)/dr) = 4:nDpo{R — e), which is precisely the leading time- 
independent term in the expression for flux in Eq. (jH)). Even though this reproduces 
correctly the leading term in the flux in Eq. (jHl), it is not entirely satisfactory since 
we are using a definition of flux which is valid only for the continuous-time diffusion 
equation and not for the discrete-time jump process. Moreover, just the knowledge of the 
stationary density profile is not enough to prove the leading time-dependent behaviour, 
namely the 1/Vi decay to the stationary flux in Eq. Thus, to rigorously prove the 
two leading terms in Eq. (jHl), we need to obtain the time-dependent solution Pn(r) of the 
discrete-time process for large n and then use this solution in an appropriate expression 
for the discrete flux (see below). This is precisely what is achieved in this paper. 
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The main achievements of this paper are: (i) to prove rigorously the result for 
the discrete flux in Eq. © for the two leading terms and then (ii) to generalize the 
results in Eqs. (0}, © and (jH)) to the case where the jump length |r — f^l from r* to 
r is not of a fixed length (as in Rayleigh flights) but is a random variable drawn from 
a distribution IV(|r — r*|) that is bounded above (with an upper cut-off < 2R), but 
is otherwise arbitrary. The stationary density profile for the general case was already 
calculated in detail in Ref. [Sj and it was shown that Eq. (jlj) is vahd for the general case 
where e has an explicit expression jH] 

"1 - fik) 



r°° dk ^ 

71 Jo k"^ 



(9) 



a^k'^/2 

where f{k) = f{x)e^^^dx is the Fourier transform of the function 

f{x) = 2n / W{u)udu , (10) 



and 0"^ = f (x) x"^ dx is the second moment of f{x) which is assumed to be finite. 
In terms of W{u), we have o"^ = (27r/3) J^^W{u) du. For the case of the Rayleigh 
flight, the jump probability W{u) = ■^;^6{u — i) in three-dimensional space translates 
to a uniform jump distribution for f{x) [TTj, 

"^'^ ^ { 0, otherwise . ^ 

Then f{k) = sin ki/{k£), and Eq. reduces to Eq. ((Tj) with e = ci. 

In addition, for the general jump distribution W{z), another remarkable universal 
result was proven jH]. It was shown that the stationary density on the surface of the 
sphere, properly rescaled by R and a, is a universal dimensionless constant [Hj 

Po.{R) R_ 1 .... 
Po ^~V2' ^ ^ 

which generalizes Eq. (jH)). In this paper, we will prove that the asymptotic expression 
for the discrete flux in Eq. (jH)) is also valid for a general jump distribution W{z) with 
the length e given by the exact formula in Eq. 0. 

Apart from the three-dimensional results mentioned above, in this paper we also 
study the one- dimensional discrete flux problem where the origin is a trap and initially 
random walkers are distributed uniformly with density po on the positive semi-infinite 
line. At every discrete time step, each walker jumps independently by a distance x 
that is chosen from a normalized probability density function f{x) with a finite second 
moment o"^ = f{x) x^ dx. Due to the recurrent nature of random walk in one 
dimension, the density approaches zero everywhere on the positive line at late times 
and the instantaneous flux to the origin also decays to zero with time. Our exact 
calculation shows that while the instantaneous flux decays algebraically as 1/ -y/n for 
large time n, the cumulative flux $™™ (i.e., the total number of particles absorbed by 
the origin up to n steps) behaves asymptotically as 



vr 12(T^ V 27rn ^ 



(13) 
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where /i4 = J^^f{x)x^dx (assuming it is finite) and e is the same quantity that 
appeared in the 3-d fiux problem in Eq. ©. The quantity e however does not appear 
in the expression for the 1-d instantaneous fiux as it does in three dimensions (at least 
up to C(n-3/2)). 

The paper is organized as follows. In section II, we set up the basic integral equation 
satisfied by the time- dependent density profile for an arbitrary jump distribution in three 
dimensions and show how one can reduce it to a one-dimensional integral equation. We 
then briefiy describe the exact solution to this integral equation first obtained in Ref. 
jH]. In section III, we use this solution to calculate $™™ and from that derive the main 
result in Eq. (jHl). In the process, we also prove analytically some other conjectures made 
in Ref. |3] on the basis of numerics. Section IV deals with the expressions for the fiux 
in one dimension. Finally we conclude in section V with a summary and outlook. Some 
of the details of the calculations are relegated to the Appendix. 



2. The Integral Equation for the density profile and its exact solution 

Let p„(r) denote the density of particles at position r outside the sphere of radius R at 
time step n. Initially we have po(^) = Po for all r > R. The density stays zero inside the 
sphere since it is a trap. A particle at position f (with r' > R) jumps in one time step 
to a new position r and the jump distance (r — f') is drawn at each step independently 
from an isotropic distribution W^(|r — f*|) that depends only on the magnitude of the 
jump length, but not on its direction. The function W{u) is normalized in the three- 
dimensional space, i.e.. 

An / W{u)u^du = 1 . (14) 
Jo 

The time evolution of the density field Pn(^) is clearly Markovian and is governed by 
the integral equation 



p„(r) = / p„_i(r-^) Wi\f- r'l) df' . (15) 

Jr'>R 

The initial condition is spherically symmetric and the evolution equation preserves this 
symmetry, hence Pn{r) = Pn{r) at all steps n. We then have 

Pn(r) = / p„-i(r') Wi\f- r^l) df' . (16) 

Jr'>R 



Using |r* — P| = + r'^ — 2r r' cos 6 where 6 is the angle between r and and 
(ir* = 2iTr'^ sin 6dr' dO in three dimensions we get 

p„(r) = 27ry dr'pn-i{r')r'^ J W (Vr2 + r'2 - 2rr'cos^) sinedO 

= — / dr' r' pn-i{r') / W{u)udu. (17) 

r JR J\r~r'\ 

Let us define a new variable Fn{r) = rpn{r)/ Pq. Then -F„(r) for r > R evolves via the 
one-dimensional integral equation 

roo rr+r' 

FJr) = 2-K dr Fn-iir') W{u)udu with Fo(r) = r . (18) 

JR J\r-r'\ 
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A further simplification can be achieved by introducing a shift, i.e., defining z = r — R 
(distance from the surface of the sphere) and writing -F„(r) = Fn{R + z) = Hn{z) where 
we have suppressed the R dependence in Hn{z) for convenience. Then Hn{z) for 2; > 
evolves via 

/•oo i-z+z'+2R 

Hn{z) = 27r dz'Hn_i{z') W{u)udu , (19) 

JO J\z-z'\ 

with Ho{z) = R + z. 

Eq. (fT^ is very general and is valid for arbitrary isotropic jump distribution 
Vr(|r — r*|). However, the reduced effective kernel in the one-dimensional integral 
equation (fTTIj) — namely f{z,z') = 2n J^l'^l,'^'^^ W (u) u du — is non-stationary since 
it depends on both r and r' and not just on its difference \r — r'\. It is difficult to 
solve this integral equation with a non-stationary kernel. To simplify further, we will 
henceforth assume that the jump distribution W{z = \r — r'l) is bounded above, i.e., 
it has a finite support with an upper cut-off -Zmax- This means W{z) = for 2; > ^max- 
If we further assume that Zmax < 2-R, it follows that one can replace the upper limit of 
integration in f{z, z') by 00 and Eq. (|T9|) simplifies to one with a stationary kernel, 

HJz)= dz' Hn^Az) fiz- z)dz' with Ho(z) = R + z, (20) 
Jo 

where f{x) is a symmetric non- negative function 

fix) = 2n / W{u)udu . (21) 



Moreover, it is easy to show from Eq. ()2H1 . via integration by parts and then using Eq. 
(HH), that f{x) is normalized to unity: J^oo fi^)^^ — 1- Thus f{x) can be regarded 
as a probability density function. For later purposes, it would further be convenient to 
break the solution of Eq. into two /^-independent parts: Hn{z) = Qn{z) + Rqn{z) 
where Qn{z) and qn{z) evolve by the same integral equation, albeit with different initial 
conditions: 

Qn{z) = / dz Qn-i{z ) f{z - z') dz with Qq[z) = z , (22) 
Jo 

Qniz) = / dz' qn-i{z') f{z - z') dz' with qo{z) = l. (23) 



Once the solutions to these two equations are known, the density profile p„(r) can be 
obtained via 

Pn{r) = ^ [Qn{r -R) + Rqn{r - R)] . (24) 

The solutions to the two integral equations in Eqs. ()22|1 and ()23|1 were obtained 
explicitly (in the Laplace domain) in Ref. jH]. We will not provide the derivation 
here, but just mention the main results. It was shown that the two following Laplace 
transforms 

/•oo / °o \ 

^,{u,s) = J^ i^Yl^qn{z)s''je-'''dz (25) 
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and 



iJu{u, s) 



have the exphcit solutions |E1 



\n=0 



i'i{u,s) = -0(0, s) 
u 



(26) 



(27) 



and 



i^ii{u, s) 



1 

u 



,s)-0'(O,s) 



[u,s) 



where the prime indicates differentiation with respect to u, 0(0, s) 



(28) 

:i - s)-'/\ and 



u,s) 



y/l — s + auJs/2 



X exp 



u 



dk 



71 Jo + k"^ 



In 



sf{k) 



(29) 



-s + sa^k^/2^ 

where as before f{k) = f{x)e^^^dx is the Fourier transform of f{x). 

Using these solutions in Eq. (j2l|l one can in principle calculate the full density 
profile at all times n. In Ref. |SI, we have shown how to obtain the stationary solution 
(the n — > oo limit) explicitly. In the next section, we show how to obtain the two leading 
terms in the expression for flux using the explicit solutions mentioned above. 



3. Calculation of the discrete Flux 

In this section our aim is to compute the expression for the discrete flux, in particular 
to prove the two leading terms in Eq. ©. To proceed, it turns out to be convenient to 
compute not the instantaneous flux ^{t) as in Eq. (jHl) but rather the cumulative flux 
^cumj^-^-j up to time t. We also replace t by the subscript n to reflect the discrete nature 
of time in this problem. From (jHl) we expect 

2{R-e)^ 



Att{R - e)Dpo 



n + 



+ 0{l) 



(30) 



where the correction term 0{1) is effectively the constant of integration. 

To prove the result in Eq. (jHUj) . we note that the cumulative flux up to time n is 
simply the total number of particles absorbed by the sphere up to time n and is given 
precisely by the "missing" mass at that time: 



IR 



47rr^[po - Pnir)]dr 



rFn{r)]dr 



Anpo / [{z + Rf -{z + R){Qn{z) + Rqn{z))]dz 



(31) 

where in line two we used Fn{r) = pn{r)/{por) and in the third line we made the shift 
z = r — R and used Fn{z + R) = Hn{z) = {Qn{z) + Rqn{z)) as explained in the previous 
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section. Here Qn{,z) and qn{z) are the solutions to Eqs. and (j2Sl) whose Laplace 
transforms are given explicitly in Eqs. and ^Tj\ respectively. So, our job is to use 
these explicit solutions in the expression in Eq. (jHlj) and obtain the two leading terms 
for large n. The rest of the section is devoted precisely to this technical task. 
Simplifying (p?T|l further, we find 

^''r = 4vrpo / {z + R)[z + R- (Q„(z) + Rqn{z))]dz 
Jo 

POO 

= 47rpo (z + R)[iz - Qniz)) + R{1 - qniz))]dz 

= 47rpo[J« + R{4'^ + Jj?^) + R'Ji'^] , (32) 

where 



= / z{z- Qn{z))dz , (33) 

JO 

J(2) = {z- Qn{z))dz , (34) 
Jo 

poo 

Jf = / z{l- qniz))dz , (35) 
Jo 

4'^ = / (1 - qn{z))dz . (36) 



One can verify that, for the case of the Rayleigh flight of jump length i the 
instantaneous flux = $™™ — $^™} calculated from the above formulas agrees with 
Eqs. (13) of Ref. W: 

$„ = ATTRDpo[a^^\n) - {e/R)b^^\n) + {R/e)a^^\n) - b^^\n)] (37) 

where 

a«(n) = 3 f\l - z)Qn{z)dz , 6«(n) = ^ /'(l - zfQ^{z)dz , 
JO I Jo 

a(2)(n) = 3 f\l - z)q^{z)dz , b'^'\n) = ^ f\l - zfq^{z)dz . (38) 
JO 2 Jo 

Now, to calculate the integrals J^^^ for a general jump- length distribution, we 

proceed as follows. We first Taylor-expand 0('U, s) defined in Eq. ()29|1 in powers of 

u for small u. This gives 



1 + a(s)M 



[a(s) + /2(s)]m + [a(s)2 + a{s)l2{s) + ^l2{sf]u^ + ... 



1 - 

(39) 



where 



0(o,.) = (l-s)-v^ (40) 

a{s) =^^=. (41) 



(1-s 






'i 


\/2(l 




1 /"a^ 




TT JO 





General flux to a trap in one and three dimensions 



9 



Note that the analytic expansion in small u in the first line of Eq. (|39p is valid up to 
0{u^) and the ... refers to non-analytic higher order terms. For our purpose, we need 
terms only up to 0{u^). 

Next, we carry out a Taylor-series expansion of -ipi and ■?/'n in powers of u: 

1 a(s) + l2(s) 



t/ji{u,s) 



u{l — s) 1 — s 

^ [a{s)' + a{s)his) + lhisY]u ^ 
1 — s 



= ^1^(1^ r^s 

^ [aisrhis)+ais)hisy + lhisr + his)]u ^ 
1 — s 

and rewrite ipi and ipn in a way that separates off the singular part in u: 

1 °° /-oo 

^iKs) = — r-E^"/ {l-qn{z))e-'^^dz (45) 

Ml — S r^r, Jo 



n=0 

oo 



1 °° /"OO 

^P,j(u, s) = - J2 s- {z-Qn{z))e-^^dz , (46) 

SO that the integrals remain convergent when m ^ 0. Expanding the exponentials e~"^ 
in powers of m, we get directly the generating functions of the desired integrals J^'^^ For 
example, s^J^^'^ follows from (minus) the coefficient of in ^l'l{u, s). Thus, from Eq. 
fjIH|l we deduce 

We extract the large-n behaviour of Jj^^^ by looking at the behaviour of the generating 
function for s — > 1. To carry out expansions to higher order, we use (see Appendix) 

^ -f:C„s", C„^2,/^+-^ asn^oo. (48) 



(1 — S)^/^ " ' V TT 4:y/7m 

We also show in the Appendix that 

Hs) - h{l) - 1^(1 - sy/' + 0{1 - .) , (49) 

and 

h{s)--^^{l-sr'/' + 0{l), (50) 

as s — > 1, where fi^ = J^^x'^f{x)dx is the fourth moment of the projected jump 
distribution /(x). 

Putting these results into ()47|1 . and using the Taylor-series expansions of (1 — sY^'^ 
and (1 — we find 
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where we have used e = — /2(1). 

Likewise, X) s^J^^ follows from the coefficient of in '^i{u, s): 



^ n m _ a{sr + a{s)l2{s) + ll2isr 



n=0 

^'s a^shjs) \h{sf 

2(1 - S)2 ^ V^(l - s)3/2 ^ 1-s ' ^ ' 

and using (jlHI), we find, 

(53, 

In the same way, s'^J^^ follows from the coefficient of in \E'ii(m, s): 



^ n ri2) _ ^i''^)Us) + \h{sf ay-Shis) , Ihjsf 



l-S v^(l-.)3/2 1_, 

which implies 



Finally, '}2s'^Jn'^ follows from the coefficient of in \E'n('U, s): 

^ n Kl) _ <S?HS) + a(^)/2(g)^ + \h{sf + /4(S) 

2^ " ~ 1 _ „ 

n=0 

ayi/2(s)^ |/2(s)^ /4(^) .... 
2(1 - sf ^ 72(1 - s)3/2 ^ 1-s ^ 1 - s ■ ^ ^ 

Note that the first term on the RHS above behaves as 

a'^sh^s) 0-25/2(1) , sy2/i4 

+ iTT-r. > (57) 



2(1 - s)2 2(1 - s)2 24a(l - sf/^ ' 
and the last term of Eq. behaves as 

74(5) S-\/2/i4 



(54) 



(5^ 



1 - s 24a(l - s)3/2 ' 
as s — s> 1, so between these two expressions, the terms containing /i4 exactly cancel out. 
With the remaining terms of (jSH)), we find 

J«~^ + ae2^ + 0(l). (59) 

The hisY term in Eq. ()56|) does not contribute to this order. 

Finally, putting the four results Eqs. dST}, (ESI), (EH) and into Eq. ((221), we 

find, after some algebra. 



27rpo(i? - e)(T^ 



n + 2(i?-e)J^ + 0(l) 



7rcr2 



(60) 



General flux to a trap in one and three dimensions 



11 



which is precisely in the form of Eq. (j3(J|) . Then $„ follows from $„ = $™™ — 
and the result is exactly in the form of Eq. © with D = a'^/2. This is the expected 
form of D given that a is the standard deviation of the one-dimensional projection of 
the three-dimensional walk. (For example, for the Rayleigh flight, this corresponds to 
the familiar formula D = [3:-) This completes the proof of our main result, Eq. 



We close this section with the following observation: If we take the difference 
between (jH^ and (jSH), we find 

tfiJ^I'-4")-^., (61) 

which implies 



exactly for all n. For the Rayleigh jump case, the rhs is ni'^/6, and in that case this is 
the cumulative equivalent to the identity a^^\n) = 1 — W^\n) observed numerically in 
Ref. Thus, we have proven that result, and shown that it can be generalized for an 
arbitrary jump distribution, with the difference given in Eq. (jU^ depending only upon 
the standard deviation a of the jump distribution. 

4. One dimension 

Let us now consider a one-dimensional system where the origin represents a trap and 
initially we have particles distributed uniformly over the positive line z > Q with density 
Po- Each particle jumps at every time step and the jump distance u is drawn from a 
symmetric jump distribution f{u). Whenever a particle crosses over the origin to the 
negative side it gets trapped by the origin. We ask similar questions as'm d = 3: What 
is the density profile Pn{z) at step n and what is the instantaneous flux at step n? 

Due to the recurrent nature of the random walks md = 1 , it is clear that eventually 
all the particles will be absorbed by the origin and hence Poo{z) = for all z > 0. This 
is a major difference from the 3-d case. In 1-d, the density profile for z > evolves in 
time via the integral equation 



Clearly, pn{z) = poqn{z) where qn{z) satisfies the integral equation (jSHI) with the initial 
condition qo{z) = 1 for all 2; > 0, and we already know its solution from section II. 
Note that g„(z) can also be interpreted [UlI as the probability that a random walker 
starting at 2; > and undergoing discrete jumps drawn from f{z) does not cross the 
origin during the first n steps, as discussed in Ref. [11]] . 

Regarding the flux, the cumulative amount is given by the missing mass at 
step n 




(62) 




(63) 




(64) 



General flux to a trap in one and three dimensions 



12 



Now, the integral J^[l—qn{z)]dz is precisely J^^ of Eq. whose asymptotic behaviour 
is given in Eq. (j^Tj) . The integral Jl^[l — qn{z)]dz also represents the expected maximum 
£'[M„] up to n steps of a random walker starting at the origin and undergoing jumps 
drawn from f{z) at every time step PU]- Thus, we have a rather nice relationship 
between the cumulative flux and the expected maximum. 



2n /i4 + 3(T 

o"W e H — 

V vr 12(73 



27m 



(65) 



where e is the same quantity in Eq. Q that appears in the expression for cumulative flux 
in three dimensions. However, while in 3-d, e retains its appearence in the expression 
for the instantaneous flux $„ = $™™ — in 1-d it disappears from at least up 

to this order. 

That is, from Eq. (jHS)), "we have 



Po (^\ 



2 L 1 



n 



-/2 _ 



n 



1) 



1/2- 



+ 



/i4 + 3a^ 
12^3 V2 



vr 



-1/2 



n 



1) 



-1/2- 



1 

2n 



P4 + 3q- 



Po- 



24^4 

-1/2 ^ 



n 



-3/2 , \ 



(66) 



/27r V 24^4 

which does not involve e. Note that the coefficient of the last term is proportional 
to the "excess," defined as 72 = f^^/cr'^ — 3. This quantity equals zero for a gaussian 
distribution. 

For the case of the uniform /(x) of Eq. (fTTj) . we have = £^/3 and fi^ = and 
from Eq. ()66|). it follows that the fiux is given by 



Po^ 




(67) 



This result proves one more formula that was conjectured in Ref. jij. For the uniform 
/(x), $n equals the poia^'^\n)/6 where a^'^\n)/6 is given in Eq. (jHHj) . and in Eq. (21c) 
of Ref. |1] is was conjectured from the numerical results that 

a(2)(n)/6 = (67r)-i/2^-i/2 ^ (i/20)(67r)-i/2^-3/2 ^ ^(^-5/2^ _ ^gg) 

Indeed, the above formula is identical to Eq. (jU7j) . Thus, we have proven Eq. (|UHj). and 
in particular verified that the coefficient in the second term is exactly 1/20. 
Note, finally, for the uniform /(x), Eq. becomes 




n 



1/2 




-1/2 



(69) 
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where c = 0.29795219 ... as given in Eq. ((7j). This result agrees with the expectation 
of the maximum of a one- dimensional walker with uniform jump density jHlCD], -E'[M„], 
which is identical to $™™ (up to the constant factor po) as mentioned above. 

5. Conclusions 

Thus, we have proven that Eq. (jHl) or its integrated form Eq. (jHn|l applies to an arbitrary 
jump distribution f{x). Here f{x) is related to the radial projection of the three- 
dimensional jump distribution iy(|r — f^l) by Eq. (|TIHl as long as W{u) has a cutoff for 
large u. The generality of the flux in Eq. (jH)) implies that there is a universal long-time 
behaviour of the intermediate particle density given by 



valid for large t and ^ r ^ V TcDt. Only the length e and diffusion coefficient D 
depend upon the details of the jump distribution. Eq. (fTOj) is equivalent to Eq. ((T)) with 
R replaced hj R — e and the complementary error function expanded to two terms. 

The integrated form of the flux of Eq. (jHUj) brings in additional constants that do 
not appear in the instantaneous form, Eq. (jHl), including the quantity e that appears 
in the one-dimensional flux, Eq. (jU^ . Indeed, e defined in Eq. © is ubiquitous as it 
appears in many terms in the expansions of the integrals J^^-* that apply to both one- 
and three-dimensional problems. 

The integrated form of the fiux also enters directly in the formula for the trap 
survival probability (the probability that an adsorbing sphere is not hit by a diffusing 
particle up to time t): S{t) = exp[— $'^"™(t)]. This quantity has been studied extensively 
and has applications to quenching of fiuorescence reactions (e.g., [121 CBl El ) • 

When the jump probability corresponds to the Rayleigh fiight, f{x) becomes the 
uniform jump distribution of Eq. ()11|) . in which case e becomes 0.29795219 . . .£ found 
previously jU For this model we have also provided proofs of the conjectured 
behaviour of the quantities a^^\n), b^^\n), a^'^\n), and b^'^\n) in Eqs. (21a)-(21d) of 
Ref. |3] and proven the conjectured identity a^^\n) = 1 — b^'^\n). We did not derive 
expressions for the higher-order corrections of Eq. (21b) and Eq. (21d) of Ref. j4j that 
involve the empirically determined constant d = 0.0270103, but these terms enter only 
in higher-order corrections to the fiux. 

We also proved the identity (jU^ for a general jump distribution. This identity 
relates the solution qn{z) with a linear initial condition to the solution Qn{z) with a 
constant initial condition, and implies a subtle relation between the corresponding one- 
and three-dimensional diffusion processes. 

In this paper we have restricted ourselves to the spatial dimensions 1 and 3. A 
natural question is: what are the expressions for the density profile Pn{f) and the 
cumulative fiux $™™ in dimensions d ^ 1,3? Is it possible to compute the Milne 
extrapolation length e for d > 37 The reason we restricted ourselves only to d = 1 and 
d = 3 can be traced back to the basic integral equation ()l(j|l . For d = 1 and d = 3, 




(70) 
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this integral equation can be reduced to a one-dimensional integral equation ((201) with 
a stationary kernel f{z — z') that is symmetric. For c? 7^ 1,3 this step cannot be done 
and finding an exact solution of the integral equation seems difficult. For example, in 
d = 2, using df* = r'dr'dO in Eq. (fTBj) one gets 

p„(r) = j'^ dr' pn-i (r') r' J^'' W ( Vr^ + r'2 - 2rr' cos O) dO. (71) 

If one performs the angular integration after the substitution = r"^ + r'^ — 2rr' cos 9 
one gets the following integral equation 



p„(r) = / rfr'p„_i(r')G(r,r') (72) 

JR 

where the kernel G{r,r') is a rather complicated function of both r and r' 

^/ r. o / r+''' dzzW(z) 

G{r,r') = 2r' / ^ ^ = . (73) 



Solving the integral equation (ff^ or the equivalent integral equations for (i > 3 remain 
a challenging problem. 



Appendix A. Proof of Eqs. (gHD and 



Proof of Eq. PHjl : For the terms that involve a{s) /(I — s), we have for s near the singular 
point at s = 1, 



(1 - s)3/2 



(1 - s)3/2 



1/2 
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1-1/2 



E(2^ + i) 

n=0 

00 

n=0 



2n\ s*^ 



n / 2 



2n 



1 ^ /2n\ ^ 



Now, for large n, 

1 /2n^ 
22" 

Therefore, we have 




1 



1 /n 1 
2n + l-- =2J- + -^ 

2 V TT 4^/7m 



+ 



(A.l) 



(A.2) 



(A.3) 



thus contributing to the n^^'^ and n~^/^ terms in the large-ra behaviour. 

Derivation of Eq. (j49j) : First differentiate Eq. (|42p (for n = 2) with respect to 
s, then put s = 1 — e where e is small. Next, inside the corresponding integrand 
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make the change of variable k = y/ey and then expand everything for small e. Use 
f{k) ~ 1 — + /i4fc'^/24. What you get as e — > is the following: 

^ ^ mi 1/2 r ^'dz 

de 127r(T3 io (1 + z^Y ' ^ ' ' 

The integral gives 7r/4. Then integrating the above with respect to e one gets 

h{e) ~ h{e = 0) + Civ^ , (A.5) 

where Ci = ^/2nJ{2Aa^). 

For J4, the integral is already singular as e — and there is no need to differentiate. 
A direct analysis of the integral gives Eq. (jSUj) . 
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